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ABSTRACT 

We have determined the three-dimensional (3D) 
structure of DNA duplex that includes tandem 
Hg"-mediated T-T base pairs (thymine-Hg"- 
thymine, T-Hg"-T) with NMR spectroscopy in 
solution. This is the first 3D structure of metallo- 
DNA (covalently metallated DNA) composed 
exclusively of 'NATURAL' bases. The T-Hg"-T 
base pairs whose chemical structure was 
determined with the ^^N NMR spectroscopy were 
well accommodated in a B-form double helix, 
mimicking normal Watson-Crick base pairs. The 
Hg atoms aligned along DNA helical axis were 
shielded from the bulk water. The complete dehy- 
dration of Hg atoms inside DNA explained the 
positive reaction entropy (AS) for the T-Hg"-T 
base pair formation. The positive AS value arises 
owing to the Hg" dehydration, which was approved 
with the 3D structure. The 3D structure explained 
extraordinary affinity of thymine towards Hg" and 
revealed arrangement of T-Hg"-T base pairs in 
metallo-DNA. 



INTRODUCTION 

The metal-mediated base pairs (the metallo base 
pairs) are currently being explored toward genetic 
code expansion (1^), development of metallo-DNAs 
(5-12), molecular magnets (13,14), electric nano-wires 
(15-19) and metal ion-sensors (20,21). Among these, 
the Hg"-sensor employing thymine-Hg"-thyniine 
(T-Hg"-T) base pair was the first successful 
application (20). 

The success of this Hg"-sensor was owing to both the 
extraordinary Ha"-thymine specificity and the thermal 
stability of T-Hg '-T base pair (22-25). The thermal sta- 
bility of the T-Hg"-T base pair was similar as those of 
normal Watson-Crick (W-C) base pairs (24). Moreover, 
the positive AS recorded for T-Hg -T base pair forma- 
tion with isothermal titration calorimetry (ITC) (24) 
indicated its pecuharity, since biomolecular complexations 
are usually hnked with negative AS values (26,27). 
However, the lack of structural data for T-Hg"-T base 
pairs in a DNA duplex prohibited rational explanation of 
this positive AS. In addition, as apparent from structure- 
based drug designs, the explanation of A 5 on structural 
basis is difficult and rarely possible. Therefore, the eluci- 
dation of entropic contributors from three-dimensional 
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(3D) structures is challenging issue in structural biology/ 
chemistry. 

The metallo-DNA with T-Hg"-T base pairs is regarded 
promising conductive nano-material. Several groups 
examined its abihty to mediate hole/electron transport, 
and weak-hole transport similar to that in normal DNA 
was actually observed (15-18). However, as mentioned 
above, the lack of structural information for T-Hg"-T 
pairs prevented rationalization of these experiments and 
tuning of conductivity in metallo-DNAs. 

The binding mode of Hg atom in T-Hg"-T base pair 
was surely determined using the '^N NMR and Raman 
spectroscopy (23,28,29). The theoretical calculations 
based on its structure suggested that the LUMOs appear- 
ing around Hg" of the T-Hg"-T base pair is distributed 
along the DNA-helical axis (29). Therefore, it is important 
to reveal mutual positioning of T-Hg"-T base pairs in the 
metallo-DNA and to confirm if the overlap of their 
LUMOs is possible or not. 

The metallophilic attraction between Hg atoms in con- 
secutive T-Hg"-T base pairs stabilizes structure of 
metallo-DNA although the Hg atoms in these metallo 
base pairs bear sizable positive charge (29-31). Only few 
observations of the metallophilic phenomenon were 
reported so far for some organometallic complexes. The 
3D structure of metallo-DNA would provide rehable basis 
for physicochemical investigation of Hg-Hg metallophilic 
attraction. 

The 3D structures of metaUo-DNAs, which are 
currently available include solely those composed of 
'ARTIFICIAL bases' (metal-chelators) (1,2,11). The 
structural information on metallo-DNAs composed of 
'NATURAL bases' is therefore very sparse. Only the 
T-Hg"-T base pair was thoroughly studied with molecu- 
lar spectroscopy (23,28,29,32-36), 3D modeling (34) and 
the crystal structure of 1 -methylthymine-Hg complex 
(37). Although the binding mode of Hg" was determined 
in these studies, 3D structure of the metallo-DNA duplex 
remained unresolved. 

The lack of structural information for the metallo-DNA 
duplex in solution made us initiate this study aiming par- 
ticularly explanation of the thermodynamic parameters 
for T-Hg"-T base pair formation. 

MATERIALS AND METHODS 

Thermal denaturation experiment 

DNA sequences used for this experiment are listed in 
Figure 1 and Supplementary Figure SI. UV spectra of 
the solutions of DNA decamers were recorded every 3°C 
for '1.2(T-Hg"-T)' and '1.3(T-A)', and every 2°C for 
'1»2(T-T)' (Figure 2). In the temperature profiles, UV ab- 
sorbances at 260 nm were plotted against temperature 
(Figure 2). The 7^ value which is dependent on the 
nearest neighbour W-C base pairs against T-Hg"-T 
base pairs was studied using a dodecamer hetero duplex: 
d(CCGCXTTVTCCG) • d(CGGAWTTYGCGG), where 
X-Y and V-W are W-C base pairs (Supplementary Figure 
SI). The effect of concentration of Hg -bound duplex 1»2 
on value was also examined (Supplementary Figure 



S2). The values were determined using a method 
described in the hterature (38). In all the thermal denatur- 
ation experiments, we confirmed that temperature profiles 
for increasing and decreasing the temperatures were iden- 
tical within the experimental error range. For further 
details, see Supplementary Material. 

NMR measurements and 3D structure determination 

NMR spectra for the ^H resonance assignments and struc- 
ture calculations were measured as described previously 
(23,39). By using the derived NMR spectra, we assigned 
all the non-exchangeable protons (39) and most of the 
exchangeable protons. The complete assignments are 
reported in Supplementary Table SI, and deposited in 
the Biological Magnetic Resonance Data Bank with acces- 
sion number 11 528. Resulting experimental constraints 
and other constraints for structure calculations are fisted 
in Supplementary Table S2. The structural constraints for 
the T-Hg"-T base pairs were generated based on the 
crystal structure of the 1 -methylthymine-Hg" (2:1) 
complex (37), i.e. N3-Hg" bond length: 2.04A and N3- 
Hg"^N3-bond angle 180°. 

Based on these structural constraints, the 3D structure 
of the DNA duplex with T-Hg"-T pairs was calculated by 
simulated annealing, using the program X-PLOR ver 3.851 
(40), based on previously reported protocols (41). From the 
calculations, 17 structures that satisfied the experimental 
constraints and covalent geometries were obtained out 
of 100 randomized structures (Supplementary Figure S3). 
Statistics for the converged structures are shown 
in Supplementary Table S2. Through the structure calcu- 
lations, the N-Hg"-N linkages of the T-Hg"-T pairs 
were maintained. This is because the pairing partners of 
each T-Hg"-T pair had already been determined in the 
same DNA sequence from the 2-bond '^N-'^N J- 
couphng across Hg" (Vnn) (23) (Figure Id). 

For further information on the NMR measurements 
and the structure calculation, see Supplementary 
Material. The structure is deposited in the Protein Data 
Bank with ID 2rt8. 

ONIOM QM/QM calculations, structural modelling and 
geometry optimization 

The structural model employed in the ONIOM QM/QM 
calculations (42); CAM-B3LYP(6-31G*, Stuttgart ECP 
for Hg):BP86(LANL2DZ) with GAUSSIAN 09 (43), 
was derived from the NMR structure of the DNA 
duplex 1»2, and is schematically depicted in Figure le 
(the G4-C17, T5-Hg"-T16 and T6-Hg"-T15 base 
pairs). The implicit water solvent was employed in all cal- 
culations. The geometry optimized structures for product 
and reactant adjusted from Equation (1) in Results and 
discussion section are depicted in Supplementary Figures 
S4 and S5. For the derivation of Equation (1), see 
Supplementary Material. In the reactant, hydrated Hg" 
bound to the DNA(T-T) while in the DNA(T-Hg"-T) 
product Hg" was completely dehydrated. The overall 
helical structure of the models was ensured by relevant 
constraints adopted from the 3D structure of Figure 3; 
only the middle base pair was geometry optimized (see 
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the legend to Supplementary Figures S4 and S5 for 
details). The A//, A5 and AG were calculated for 
r= 298.15 K and standard pressure within the rigid- 
rotor harmonic-oscillator approximation, 5 was 
composed of translation, rotation and vibration contribu- 
tions. For ONIOM QM/QM calculations, see also 
Supplementary Material. 

Preliminary X-ray diffraction analysis of the DNA duplex 
and modelling of the DNA duplex 1»2 with an 
experimental Hg-Hg distance constraint 

To obtain experimental Hg"-Hg" distance in a DNA 
duplex, a preliminary X-ray diffraction analysis was per- 
formed. For this purpose, the DNA dodecamer was co- 
crystallized with Hg(C104)2 by the hanging-drop vapour 
diffusion method at 4°C. Preliminary X-ray data collec- 
tions were performed with synchrotron radiation 
(1 = 0.98 A) at BL17A in the Photon Factory (Tsukuba, 
Japan). Two Hg" atoms were found at coordinates (x, j, 
z) = (1.1, 0.0, 7.2), (1.6, 0.9, 10.3) (Supplementary Figure 
S6) using the heavy-atom-search procedure of the 
program AutoSol from the Phenix suite (44^6). Then, 
the Hg"-Hg" distance was determined as 3. 3 A which is 
also consistent with its theoretical values, 3.28-3.52A (30). 
Based on these facts, the model structure of the DNA 
duplex with T-Hg"-T pairs was also calculated by using 
rigid body minimizations and following normal energy 
minimizations under the Hg"-Hg" distance constraint 
(3. 3 A). The derived model structure is shown in 
Supplementary Figure S7. For details on the preliminary 
X-ray diffraction analysis and the modelhng studies, see 
Supplementary Material. 

RESULTS AND DISCUSSION 

Thermal denaturation experiment 

Two non-self-complementary DNA duplexes presented in 
Figure 1 were chemically synthesized. The thermal 
stabilities (r^) of duplex 1»2 with T-T mismatches and 
T-Hg"-T base pairs were evaluated and compared with 
that of reference duplex 1»3 containing W-C T-A base 
pairs (Figure 2). 

The I'm value for duplex 1»2 increased from 25°C to 
54° C upon adding Hg' . Interestingly, the DNA duplex 
1»2 with T-Hg"-T pairs was more stable than its refer- 
ence duplex 1»3 with W-C base pairs. The T-Hg"-T base 
pair is therefore more stable than the W-C base pair in the 
sequence context of the DNA duplex 1»2. We also studied 
the effects of T-Hg"-T nearest neighbour base pairs on 
value (Supplementary Figure SI) and the highest sta- 
bihty was observed for the closely related sequence with 
the Hg"-bound DNA duplex 1»2. 

Structure determination 

The DNA duplex 1»2 was selected for structure determin- 
ation because: (i) the chemical structure of T-Hg"-T was 
determined with Vnn^ 2-bond '^N-'^N /-couphng across 
Hg"-mediated linkage (Figure Id) for the same sequence 
of DNA oligomer (23), and (h) the closely related 
sequence was thermally most stable (Supplementary 
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Figure 1. The DNA sequences and the T-Hg"-T base pair, (a) The 
DNA oligomers 1, 2 and 3. (B) The DNA duplex 1»2 with residue 
numbers. The definition of left and right segments is depicted, 
(c) The control DNA duplex 1»3 with residue numbers, (d) The 
reaction scheme for T-Hg"-T base pair formation (the proton-Hg" 
exchange reaction) and 2-bond '""N-'^N /-coupling (Vnn) (23). (e) 
The schematic representation of the model used in ONIOM QM/QM 
calculation of thermodynamic parameters. The DNA(T-T) and 
DNA(T-Hg"-T) stand for Hg"-free and Hg"-bound three base- 
paired (3 bp) duplex. 
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Figure 2. The temperature profiles of UV absorbance at 260 nm. The 
vertical axis is a relative absorbance normalized between absorbances at 
the lowest and the highest temperatures. Blue diamonds: the DNA 
duplex 1»2 in the absence of Hg". Pink squares: the DNA duplex 
1»2 in the presence of Hg". Green triangles: the DNA duplex 1»3. 
The values for these profiles are given by the red characters with 
arrows. 



Figure SI). We then recorded NOESY spectra of Hg"- 
bound DNA duplex 1»2 (39), and generated distance con- 
straints from NOESY spectra published in the reference 
(39) (Supplementary Table S2). 

In total, 17 structures that satisfied the NOE constraints 
were obtained (Supplementary Figure S3 and Table S2). 
All the derived structures were normal B-form duplexes 
(Suppleinentary Figure S3). The 3D structure of the 
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(c) 




G11 



Figure 3. The 3D structure of Hg"-bound DNA duplex 1»2. (a) The 
side view perpendicular to helical axis, (b) The top view along helical 
axis, (c) The space-filling model of the middle 3-bp DNA segment 
including G4-C17, T5-Hg"-T16 and T6-Hg"-T15 pairs sketched out 
in Figure le. The Hg atoms are depicted always as red balls. The Hg- 
Hg distance derived solely with NOEs was ~4A. When we applied the 
Hg-Hg distance constraint at 3.3 A reflecting our X-ray diffraction 
analysis of a DNA duplex with tandem T-Hg"-T base pairs (53) 
(Supplementary Figure S6 and Suppleinentary Material), the derived 
model structure of duplex 1»2 (Supplementary Figure S7) was consist- 
ent with the NOE constraints. 



duplex with the lowest energy is shown in Figure 3. The 
T-Hg"-T pairs are well stacked without distorting the 
duplex, which indicates that they are accommodated 
into the DNA duplex in a similar manner to canonical 
W-C base pairs. 

In addition, due to the local topology of the T-Hg"-T 
pairs, their C 1 '-C 1 ' distances are shorter by ~ 1 A than those 
found in W-C base pairs. Nonetheless, the difference was 
within the structural variation of B-form DNA duplexes. 
Therefore, the T-Hg"-T pairs structurally mimic W-C 
base pairs without any significant perturbation of the 
double helix. This fact most likely explains why the DNA 
polymerase can incorporate thymine (T) against T in the 
template strand via the formation of T-Hg"-T pair (47,48). 

Closer look along the helical axis of DNA duplex 1»2 
revealed perfectly ahgned Hg atoms (Figure 3b). The 
space-filling model of the respective part of DNA duplex 
1»2 further revealed that Hg atoms are shielded from bulk 
water (Figure 3c). The weU-stacked metaUo base pairs 



with the Hg"-Hg" distance at 4.03^. 17A and narrow 
04— 04/02— 02 spacing exclude any possibihty for bulk 
water to penetrate into proximity of Hg atoms. The rela- 
tionship between 3D structure and the thermodynamic 
parameters will be discussed later. 

Based on the 3D structure, we confirmed the theoretical 
prediction by Voityuk (49) and us (29) that overlap of 
their LUMOs of Hg atoms in the metallo-DNA is 
possible. This implies that the metallo-DNA duplex 
could be effective route for (an) excess electron(s). 
Furthermore, the well-stacked arrangement of T-Hg"-T 
pairs suggests that interaction of Hg atoms is not repul- 
sive, which supports existence of the Hg-Hg metallophihc 
attraction inside metallo-DNA. Recently, existence of the 
metallophilhc attraction between heavy metals in metallo- 
DNAs was theoretically proposed for the Ag'-mediated 
imidazole-imidazole base pairs (50) and the T-Hg"-T 
pairs (30,31), and the 3D structure presented here is an 
additional indicative of such a newly proposed attractive 
force between heavy metals. 

The relationship between 3D structure and the 
thermodynamic parameters 

Next, we considered relationship between the 3D structure 
and the thermodynamic parameters for the T-Hg"-T 
formation determined by Torigoe's group (24,25), 
which showed positive A5 and negative A// (Table 1). 
Based on the 3D structure, the reaction for T-Hg"-T 
base pair formation can be written as follows. See 
also Supplementary Material for detailed derivation of 
Equation (1). 



DNA(T - T) + [Hg(H20)„(0H)]+ 

DNA(T-Hg"-T)+H30++nH20 



(1) 



where the DNA(T-T) and DNA(T-Hg"-T) stand for 
DNA duplex with T-T mismatch and that with T-Hg"- 
T pair, respectively (Figure le). In the Equation (1), we 
considered (i) the imino proton (H^)-Hg" exchange 
reaction upon T-Hg"-T base pair formation 
(Figure Id); (ii) the dehydration of Hg" cation during 
reaction; and (in) the ^K^ = 3.4 for Hg"-aqua complex 
(51) that implies existence of hydroxy-ligand of Hg". 

Upon the Hg"-binding to T-T mismatch, number of 
water molecules initially coordinated to Hg" were 
released to bulk. Accordingly, the dehydration of Hg" 
should yield the entropy increase following the thermo- 
dynamic assumptions. Such positive is known as 
dehydration entropy. However, the complete dehydration 



Table 1. Experimental and theoretical thermodynamic parameters 





Aiy'/kcal/mol 


A57cal/mol/K 


AG°/kcal/mol 


Reference 


Experimental (ITC)" 


-3.85 ± 0.18 


13.1 ± 0.65 


-7.76 ± 0.19 


(24) 




-4.76 ± 0.13 


10.6 ± 0.84 


-7.91 ± 0.12 


(24) 


Theoretical'' 


-4.04 


14.2 


-8.27 


This work 



"In reference (25), thermodynamic parameters possessed much larger standard deviations. Therefore, only the 
precise data from reference (24) were shown in table. ''Calculated values are based on Equation (1) (see 
Supplementary Figures S4 and S5, and Supplementary Methods). AG° values are given at 298. 1 5 K. 
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in this case has been only rarely vahdated experimentally. 
In summary, based on the 3D structure, one contributor 
to the AS was identified as dehydration entropy owing to 
the complete dehydration of Hg". 

The 3D structure of Hg"-bound DNA duplex 1»2 
enabled calculation of the A^* and AH with the ONIOM 
QM/QM method. Using the 3-bp DNA segment 
(Figures le and 3c) derived from the 3D structure of 
DNA duplex 1»2, the models of product (Supplementary 
Figure S4) and reactant (Supplementary Figure S5) were 
constructed foUowing the Equation (1). Based on these 
structures, the thermodynamic parameters were calculated 
(Table 1). The calculated AH (-4.04kcal/mol) and A^* 
(14.2cal/mol/K) agreed with the experiment (Table 1). 
From the result, not only the absolute values of AH and 
A^, but also the positive sign for was reproduced by 
theory. In addition, the calculated thermodynamic param- 
eters in this work were consistent with those derived pre- 
viously with different protocols for the complete reaction 
pathway describing formation of T-Hg"-T base pair (52). 
The determined 3D structure rationally explained the 
thermodynamic parameters. 

CONCLUSION 

The first 3D structure of metallo-DNA composed 
exclusively of 'NATURAL' bases and containing 
tandem T-Hg"-T base pairs was determined in solution. 
The positive AS recorded for T-Hg"-T base pair forma- 
tion was explained as Hg"-dehydration entropy on the 
structural basis. The 3D structure rationally explained 
the specific Hg" affinity toward T-T mismatch and 
unveiled the 3D arrangement of the metallo base pairs 
in the DNA duplex. 
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